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ABSTRACT 

Based on a three-dimensional model of an early star-forming galaxy, we explore 
the evolution of the sub-millimeter brightness. The model galaxy is employed from 
an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori 
& Umemura, where the star formation rate (SFR) is ~ 10 Mq yr~^ at tagc ^ 0.3 
Gyr and several Mq yr~^ at tago > 0.3 Gyr. The former phase well reproduces 
the observed properties of Lyman alpha emitters (LAEs) and the latter does Lyman 
break galaxies (LBGs). We solve the three-dimensional radiative transfer in the clumpy 
interstellar media in this model galaxy, taking the size distributions of dust grains into 
account, and calculate the dust temperature as a function of galactic evolutionary 
time. We find that the dumpiness of interstellar media plays an important role for 
the sub-millimeter brightness. In the LAE phase, dust grains are concentrated on 
clumpy star-forming regions that are distributed all over the galaxy, and the grains 
can effectively absorb UV radiation from stars. As a result, the dust is heated up to 
Tdust ^ 35 K. In the LBG phase, the continuous supernovae drive dust grains far 
away from star- forming regions. Then, the grains cannot absorb much radiation from 
stars, and becomes into a cold state close to the CMB temperature. Consequently, 
the dust temperature decreases with the evolutionary time, where the mass-weighted 
mean temperature is Tdust = 26 K at iage = 0.1 Gyr and Tdust = 21 K at tage = 
1.0 Gyr. By this analysis, it turns out that the sub-millimeter brightness is higher in 
the LAE phase than that in the LBG phase, although the dust-to-gas ratio increases 
monotonically as a function of time. We derive the spectral energy distributions by 
placing the model galaxy at a given redshift. The peak flux at 850 ^m is found to be 
Ss50 ~ 0.2 — 0.9 mJy if the model galaxy is placed at 6 ^ z ^ 2. This means that 
ALMA can detect an early star-forming galaxy with SFR of ~ 10 Mq yr~^ by less 
than one hour integration with 16 antennas. 

Key words: radiative transfer - ISM: dust, extinction - galaxies: evolution - galaxies: 
formation - galaxies: high-redshift - methods: numerical 



1 INTRODUCTION 

Exploring the properties of early star-forming galaxies is one 
of key themes to elucidate the galaxy formation process. The 
fashion of star formation in galaxies can significantly change 
dust properties. In young galaxies, type-II supernovae ere- 
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In addition, UV radiation from young stars is effectively ab- 
sorbed by dust and therefore determines the temperature 
of dust grains. Besides, the dust can affect the star forma- 
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tion efficiency itself via hydrogen molecule formation and the 
cooling by thermal emission. Accordingly, dust properties 
such as amount, temperature, size and composition should 
tightly correlate with the star formation history. Hence, un- 
derstanding the dust properties is a significant issue in the 
study of galaxy evolution. 

So far, the sub-millimeter (sub-mm) flux from dis- 
tant galaxies has been detected at redshifts of 0.3 — 
4. Th ese galaxies are calle d sub- m m galaxies ( S MGs) 
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ited to extremely luminous starburst galaxies, where the star 
formation rate (SFR) is ~ 1000 M0 yr~^, fro m the restric- 
tion by the detection limit of the facility (e.g.. [Borvs et all 
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dust properties in normal star-forming galaxies (SFR < 
100 Mq yr~^) at high redshifts (which mean z > 3 in this 
paper) have not been hitherto understood well. 



Recently, some observations are dedicated to detect 
the sub-mm flux from early star-forming galaxies (SFR 
~ 0.1 - 100 Mq yr~^) like Lyman alpha emitters (LAEs) 
and Lyman bre ak galaxies (LBGs) by using newly deve l- 
oped faci lities (iTamura et al. l2009l : IStanwav et al.l l2010t ). 
Although iTamura et al. ( 200£ ) did not detect notable sub- 
mm flux from LAEs and LBGs in the SSA22 region within 
the detection limit of Shoo ^ 2 mjy , they discovered sub - 
mm sources in the vicinity of LAEs. IStanwav et al.l (|2010l ) 
showed that the upper-limit of sub-mm flux for LBGs at 
z ~ 5 is SgTo < 0.85 mJy from stacked data. Hence, with 
recent sub-mm telescopes (e.g. Hershel, AzTEK), it seems 
very difflcult to detect sub-mm flux from high-z star-forming 
galaxies with SFR< 100 Mq yr~^. Moreover, the determi- 
nation of dust temperature is a hard task in the obser- 
vational study, because it requires the flux data at many 
wavelengths in sub-mm band. Hence, one often assumes the 
dust temperature to speculate spectral energy distributions 
(SED). Then, the estimation of Lm, by integrating specu- 
lated SED suffer from the uncertainty of the assumed dust 
temperature. The uncertainty leads to the misestimation of 
SFR and dust amount. A theoretical model may help to 
relate the observed sub-mm flux with the dust properties 



in star-forming galax ies. Such attempts have been made by 
several authors fe.g.. Fliikagi. Vansevicius. fc Arimoto l2003l : 



iDaval. Hirashita. fc FerraraH201Gl '). lTakagi et all (|200a ') con- 
structed an SED model of dusty star-forming galaxies 
by solving t he ra diative transfer in spherical symmetry. 
iDaval et aLl(|201Gt ). based on cosmological SPH simulations, 
employed homogeneous dist ributions of stars and d ust in 
an identified LAE. Recently, iMori fc Umemural (|2006l ) have 
shown through ultra-high-resolution numerical simulations 
that interstellar media in early star-forming galaxies become 
quite clumpy as a result of multiple supernovae. Hence, 
to construct the SED for early star-forming galaxies, it is 
imperative to solve three-dimensional radiative transfer in 
clumpy interstellar media. 



Here, we theoretically explore the dust amount and 
temperature in star-forming galaxies by solving three- 
dimensional radiative transfer, based o n a chemodynamic 
simula tion of a primordial galaxy by iMori fc Umemural 
l)2006h . Since the model galaxy well reproduces the observed 
properties of Lyman alpha emitters (LAEs) at tage ;$ 0.3 
Gyr and Lyman break galaxies (LBGs) at tage > 0.3 Gyr, 
we can predict the sub-mm brightness in the LAE and LBG 
phases in early star-forming galaxies. In the present analy- 
ses, the cosmological parameters are assumed to be Hq — 70 
km s"^ Mpc~\ Qm = 0.3 and = 0.7. In §2, the model 
and numerical method are described. In §3, the results on 
dust properties in early star-forming galaxies are presented. 
In §4, the detectability by ALMA is discussed. §5 is devoted 
to the summary. 



2 MODEL & METHOD 
2.1 Model galaxy 

Our model galaxy is a supernova-dominated star-forming 
galaxy simulated by an ultra-high-resolution (1024^ fixed 
Cartesian grids) chemodynamics calculation coupled with 
the collisionless dynamics of star particles. The simulation 
pursues the early evolution (< 2 x 10® years) of a primeval 
galaxy as an assemblage of sub-galactic condensations of 
5.0 X 10® Mq, building up a system wi th the total mass of 
10" Mq [see iMori fc Umemural (120061 ) for the simulation 
details] . 

The star formation rate (SFR) of the model galaxy is 
11 Mq yr"^ at fage = 0.1 Gyr, 10 Mq yr"^ at tagc = 0.3 
Gyr, 8 Mq yr"^ at tagc = 0.5 Gyr, and 5 Mq yr~^ at 
iage ~ 1 Gyr. The metal and dust in the galaxy is enriched 
by the continuous type II supernovae due to vigorous star 
formation, and exhi bits complex inhom ogeneous dust distri- 
butions (figure 1 in lYaiima et al.ll2009l ). 

According to the Lyman a luminosity, the early evolu- 
tionary stage is divided into two phases: one is the LAE- 
phase (fago ;$ 0.3 Gyr) and the other is the LBG-phase 
(tago > 0.3 Gyr). Most of Lyman a emission comes from 
the cooling radiation by interstellar gas, and the luminosity 
reaches 2.0 x lO"*^ erg s"^ at tage = 0.1 Gyr and 1.6 x lO"*^ 
erg s~^ at tago = 0.3 Gyr, respectively. They nicely match 
the observed luminosity of LAEs and also well resemble 
LAEs with respect to other properties. At fago ^ 0.5 Gyr, 
the Lyman a luminosity quickly declines to several 10**^ 
erg s~^ that is lower than the detection limit, and the 
SED is dominated by stellar continuum. This phase appears 
like LBGs. The model galaxy at this phase has stellar age 
(~ 2 - 6 X 10** yr) and mass (~ 6 - 9 x 10® Mq) similar 
to typical LBGs at z ~ 3, while it is older and more mas- 
sive than typical LBGs at z > 5 (|Verma et al.ll2007l '). On 
the other hand, a part of LBGs at z > 5 have been de- 
tected by Spitzer IRA Q, and show the s i milar age and mas s 
to our model galaxy (|Evles et all l2007l : IStark et al] 120091 ). 
Therefore, our model galaxy at a later phase is probably 
corresponding to typical LBGs for z ^ 3, and a massive 
subset of LBGs for z > 5. 



2.2 Radiative transfer 

The radiation from young stars propagates in the highly 
inhomogeneous interstellar media containing dust. We com- 
pute the three-dimensional radiative transfer (RT) of stellar 
radiation to derive the dust temperature. The RT calcula- 
tions are done as the post-processing for each evolutionary 
stage of model galaxy. For the purpose, the data of the hy- 
drodynamic simulations are coarse-grained into 128^ Carte- 
sian grids. 

The RT scheme used in this paper is the Authentic Ra- 
diation Transfer (ART) method which is ori ginally devel- 
oped by iNakamoto. Umemura. fc Susaj (|200ll). The proce - 
dure is basically the same a s that in Yaiima et al.l (|2009l ): 
lYaiima. Choi, fc Nagamine I (|201lf ). Our ART method is 
based on the long-characteristic method. Usually, the short- 
characteristic method is computationally less costly than 
the long-characteristic method, by an order of A'^ which is 
a grid number in the linear dimension. However, the short- 
characteristic method suffers from numerical diffusion effect. 
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The ART method is devised to reduce the computational 
cost to a level similar to the short-characteristic method with 
keeping the accuracy equivalent to the long-characteristic 
method. Hence, the present method allows us to calculate 
the transfer of radiation from a large number of sources. 

In this work, the RT equation is solved along 16384 rays 
with uniform angular resolution from each star particle. The 
number of star particles is 3 — 8.5 x 10*. Hence, we carry 
out the RT calculation for ~ 10*~^ rays for each snapshot. 
We map individual star particles to nearest grids, and set 
radiation rays in an isotropic fashion from each star particle. 



2.3 Dust model 

We distribute the interstellar dust in proportion to the 
metall icity, assuming the size distributio n of dud/dad oc 
a^^-^ (|Mathis. Rumpl. fc Nordsiecklll977l '). where Od is the 
radius of a dust grain. We suppose the grain size in the 
range of 0.1 — 1.0 /im as our fiducial model. The dust mass 
in a cell is calculated by a following simp le relation between 
metallicity and dust (|Draine et al.ll2007l ). 



TTld 



0.01m,, 



(1) 



where rud, rUg, and Z are the dust mass, gas mass, and 
metallicity in a cell. The density in a dust grain is assumed 
to be that of silicate-like grains, 3gcm~^. The optical depth 
by dust is given by 



drdust = / Q{ad,iy)TTal ( 



da^ds, 
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where Q{ad,v) and rid are the absorption coefficient factor 
and the grain number density, res pectively. We adopt the 
Q- value of silicate grain derived in iDraine fc Led (|l984 ) . 

Of course, there are some options for the dust model. 
The pre sent dust model is m otivated by the supernova dust 
model in lNozawa et al.l (|2007l l. In their model, the small dust 
grains of < 0.01 pm are readily destroyed in shock wave, so 
that the typical si z e bec omes ~ 0.1 /xm. On the other hand, 
iTodini fc Ferraral (|200lh pointed out that the typical size 
could be reduced to ~ 0.01 /xm for the first grains in the 
early universe. The number density of dust grains increases 
with decreasing dust size for the constant amount of dust, 
while the absorption cross section decreases with the dust 
size. In Figure [T] the extinction per unit dust mass is com- 
pared between ~ 0.14 ^m and ~ 0.01 ^m dust. As seen in 
Figure [T] the absorption efficiency does not change appre- 
ciably with changing the dust size. Therefore, the infrared 
luminosity and dust temperature are not so sensitive to the 
dust size. Actually, in a test calculation for 0.01 /im dust, 
the temperature decreases by just a few K. The diff'erence in 
the peak fiux at 850 /xm is a few per cent, because the peak 
of thermal emission shifts toward 850 ^m with decreasing 
temperature. 

On the other hand, the Calzetti's law 
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0.657(-2.156 + 1.509/A - O.IQS/A^ + O.Oll/A'') + 1 
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(|Calzetti et al.llioool ') is frequently used for dust extinc- 
tion in local galaxies. In Figure[2l the extinction curve of our 
dust model is compared to the Calzetti's law. Our model 
is somewhat fiatter than the Calzetti's law in UV-optical 
range. However, even if we use Calzetti's extinction curve, 
the infrared luminosity does not significantly change from 
our model. Actually, by a test calculation, we find that the 
relative difference is twenty per cent for IR luminosity, and 
the difference in dust temperature is ~ a few K. 

2.4 Radiative equilibrium 

We evaluate the dust temperature Tdust by solving the ra- 
diative equilibrium between heating by photo-absorption (F) 
and cooling by thermal emission (A), which are given by 



r = 
A = 
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where and fccMB are the mean intensity of stellar radia- 
tion and the heating term by CMB radiation, respectively. 
We simplify the right side of this equation as follows (e.g., 
lEvanslll994l ), 



// 



'i-K'^alB^{Tdust)Q{a.d,i^)^^r^ da dv ~ 4n'^alndQ(rTduBt, 
dad 



(5) 

where a is the Stefan-Boltzmann constant, ad is the mean 
dust size weighted by the size distribution function, 



ad 



(6) 



The Q is a mean Q-value weighted by the Planck function 
/ Bd(Tdust,u)Q{dd,i')di' 



Q (Tdust) — r r> rrri \J 

J Bd{Tduat,l')dl' 

Then, we obtain the dust temperature by 

1/4 

Td 
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3 RESULTS 

3.1 Evolution of dust component 

In Figure [3l we present the time variation of dust compo- 
nent as a function of the galactic evolution time. Figure[21[a) 
shows the dust-to-gas mass ratio (D). The ratio increases 
monotonically with the galaxy evolution owing to continu- 
ous supernovae. D is 4 x 10~* (corresponding to the metal- 
licity Z = 4 X 10~ Zq, where Zq is the solar metallicity) at 
0.1 Gyr and reaches 0.7 x 10"^ {Z = 0.7^©) at 1.0 Gyr. 

Figure [3{b) shows the total dust mass in the model 
galaxy. The dust mass is 1.6 x 10 Mq at 0.1 Gyr, and then 
(6 — 7) x 10^ Mq at ^ 0.3 Gyr. In contrast to the monotonic 
increase of metallicity, the total dust mass does not change 
greatly at > 0.3 Gyr. This trend comes from the fact that 
the continuous energy input by multiple supernovae results 
in the blowout of dusty gas from the halo. The fraction of 
escaped gas to the initial amount becomes ~ 50 per cent at 
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3.2 Dust temperature 
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Figure 1. Dust-absorption cross-section per unit dust mass. 
Solid and dash lines sliow tlie cross section of 0.14 /jm and 0.01 
silica te dust, respe c tively . The Q-value of silicate dust is taken 
from lDraine fc Led lll984l) . 







Our model 
Colzetti's law 



0.2 



0.£ 



1.0 



0.4 0.6 
A (jtim) 

Figure 2. The dust extinction curve. Solid and d ash lines show 
the d ust model in this paper and Calzetti's law ijCalzetti et al.l 
l2000t) . respectively. 



Our model shows that the metaUicity reaches the level 
of Q.2Z(^ < Z < OAZpf in th e late phase of LAE. Very 
recently, iNakaiima et al.l (|201ll ') assessed the metaUicity of 
LAEs to be Z > Q.16Zq, which is significantly higher than 
was previously thought for LAEs. This is concordant with 
our results. 

Figure shows the radiation energy absorbed by 

dust, and Figure Efd) shows the absorbed fraction to the 
total radiation energy emitted by stars. The model galaxy 
has a blue SED owing to young stars, and the absorption 
efficiency by dust decreases steeply with increasing wave- 
length. Thus, most of energy absorbed by dust is attributed 
to UV-optical continuum radiation from stars in the range 
of A ~ 1000 - 6000 A, and strongly depends on SFR and 
dust amount. Figure[3l^c) shows that the absorbed energy by 
dust is peaked at 0.3 Gyr, although the metaUicity increases 
monotonically. This is because the dust is distributed com- 
pactly around star forming regions in an early phase, and 
hence the UV photons from stars are effectively absorbed. 



Figure |4] shows the dust temperature distributions, when the 
model galaxy is placed at z = 3. The dust temperature is 
Tdust ~ 35 K near star-forming regions, while Tdust ~ 10 K 



far from there. In an early phase (fa 



0.1 - 0.3 Gyr) 



the high-temperature regions are distributed extensively ac- 
cording to the distributions of star-forming regions. In a 
later phase (fago = 0.5 - 1.0 Gyr), the distributions of high- 
temperature dust are confined into the central regions as a 
result of the dynamical evolution of star-forming regions. 

We statistically study the dust temperature distribu- 
tions. Figure [5] shows the mass fraction of dust in a given 
temperature range. In an early stage, dust is confined in each 
sub-halo, and therefore it is distributed near star-forming re- 
gions. Thus, at fagc =0.1 Gyr, a part of dust has high tem- 
perature of Tdust ^ 35 K, and the mass fraction is broadly 
distributed in the range of Tdust = 10 — 40 K. Thereafter, 
the dusty gas is blown away by supernova explosions and 
is distributed diffusely in the extended regions. Then, a 
part of dust, which is far from young stars, cannot absorb 
much radiation energy, and therefore becomes in a cooler 
state. Hence, in a later phase, dust with high temperature 
of Tdust ^ 35 K disappears, and the dust temperature falls 
in the range of Tdust ~ 10 — 30 K. 

Figure [H] shows the mean temperature as a function of 
evolution time. Filled and open circles are the mass weighted 
mean temperatures and the energy weighted ones, respec- 
tively. The mass weighted mean is ~ 24 K in the LAE phase, 
while ~ 20 K in the LBG phase. The energy weighted mean 
is ~ 31 K in the LAE phase, while ~ 26 K in the LBG 
phase. The high-temperature dust emits thermal radiation 
more effectively than colder one owing to the higher Q- 
value and the strong temperature-dependence of emissivity. 
Hence, the energy weighted mean temperatures are higher 
by ~ 5K than the mass weighted ones. The mean temper- 
atures decline with the evolution time especially in early 
phases. This is intimat ely relevant to the su b-mm brightness 
of the galaxy. Recently. iHwang et al.l (|2010l ) have studied the 
dust properties of galaxies with Lir = 10^^ Lq in the red- 
shift range 0.1 < z < 2.8, and found the dus t temperature to 
be Tdust ~ 20 - 50 K. lAmblard et Zl (|2010l ) have found that 
the average dust temperature is 28±8 K for submm galaxies 
with the average redshift of 2.2 ± 0.6. These temperatures 
are comparable to the dust temperature in our calculations. 



3.3 Spectral Energy Distributions 

We compute the intrinsic spectral energy distributions 
(SED) of the stellar component in the model g alaxy by using 
the p opulation synthesis code PEGASE v2.0 (|Fioc fc Roccal 
Il997h . We place the model galaxy at several redsh ifts. The 
initial mass function (IMF) is assumed to be that bv lSalpeterl 
(|l955h in the mass range of 0.1 — 50Mq. Also, the effects by 
the age and metaUicity of stellar population are incorporated 
by interpolating the table generated by PEGASE. The dust 
temerature is primarily determined by the UV continuum 
from OB stars with > 2Mq. For the Salpeter IMF in the 
range of O.l — 50Mq , the mass fraction of stars with > 2Mq 
is 0.267. If we suppose the range of O.l — WOMq, the frac- 
tion is 0.287. Hence, the upper bound of mass range is not 
so significant on determining the dust temperature. 
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Figure 3. upper left : Dust-to-gas mass ratio as a function of evolution time, upper right : Total dust mass as a function of evolution 
time. The dust mass is estimated in proportion to metallicity. lower left : Radiation energy absorbed by dust as a function of evolution 
time (filled circles). Open circles are intrinsically radiated energy from stars, lower right : The absorbed fraction to the total radiated 
energy as a function of the dust mass. 



In this work, at first, we derive tlie intrinsic SED of stel- 
lar component by PEGASE, and correct the SED by incor- 
porating the dust extinction through three-dimensional RT 
simulations. Then, we evaluate the dust temperature and 
the thermal emission following the equation (7). Finally, we 
get the emergent SED by combining the dust absorption- 
corrected stellar SED with the thermal emission from dust. 
Fig. [7| shows the resultant SED of the model galaxy in 
the observed frame after making the K-correction. As well 
known, A > 1 mm, the flux density does not decrease with 
redshift owing to the effect of negative K-correction. 



3.4 Sub-millimeter brightness 

Here, we see the time variations of the sub-millimeter bright- 
ness, focusing on the flux at 850 nm. Figure [8] shows the 
total flux at 850 /im in observed frame (Sgso)- Interestingly, 
the sub-millimeter brightness is peaked at 0.3 Gyr (LAE 
phase), independent of the assumed redshift. As already 
shown in Figure [JJ a), the dust-to-gas ratio (metallicity) in- 
creases monotonically with time. But, the sub- millimeter 
brightness does not increases according to the increase of 
metallicity. This is because, in the LAE phase, dust grains 
are concentrated on clumpy star-forming regions and there- 
fore the dust grains can be effectively heated up by stel- 
lar radiation. On the other hand, the continuous supernovae 
drive dust grains far away from star-forming regions at 0.5 
Gyr (LBG phase). Then the grains cannot absorb much ra- 



diation from stars and turn into a cold state. As a result, 
the sub-millimeter brightness declines. 

The relative sub-millimeter flux at each epoch depends 
on the redshift z. If the galaxy is at 2 < 3, the total flux 
at 0.1 Gyr is smaller than that at 0.5 Gyr and 1.0 Gyr, 
although the total absorbed energy is nearly equal. This 
is due to the fact that the flux of 850 /xm is sensitive to 
the dust temperature. The spectral shape of dust emission 
is basically that of the black body, although the Q-value 
affects the shape to some degree. The temperature, which 
exhibits a peak at 850 fim in the observed frame, is 13.65 K 
for a galaxy at 2; = 3. As seen in Figure [S] the mass- weighted 
mean temperature declines with the evolution time, and the 
temperature at 0.5 Gyr and 1.0 Gyr is closer to 13.65 K 
than that at 0.1 Gyr. Therefore, we conclude that the low 
sub-mm flux in an earlier LAE phase is not attributed to the 
low dust amount but rather to the higher dust temperature. 
For higher galaxy redshifts, the absolute flux decreases of 
course with increasing redshift. But, the difference of the 
flux between 0.1 and 0.5 (or 1.0) Gyr becomes smaller, since 
the 850 fim in the observed frame is corresponding to shorter 
wavelengths. 



4 DETECTABILITY BY ALMA 

Here, we assess the feasibility to detect early star-forming 
galaxies with SFR of ~ 10 Mq yr~^ by Atacama Large Mil- 
limeter /submillimeter Array (ALMA). In the present sim- 
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Figure 4. The distributions of dust temperature wiien tiie model galaxy is placed at 2: = 3. The color shows dust temperature. 
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Figure 9. The map of brightness temperature at A = 850 ^m in the model galaxy, assuming the galaxy to be located at 2 = 3. The 
linear scale of each panel is 16.6". Upper panels show the map with the original spatial resolution which corresponds to the angular 
resolution of 0.13". Lower panels are the map coarse-grained with the angular resolution of 0.65". 
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Figure 7. Spectral energy distributions of the model galaxy in 
the observed frame after the K-correction is made. Different line 
style corresponds to a different redshift at which the model galaxy 
is placed (solid: z = 2, dot: z = 3, dash: z = 4, long-dash: z = 5 
and dot-dash: z = 6). 
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Figure 6. The mean temperature of dust as a function of evo- 
lution time. Open and filled circles are the energy weighted mean 
and the mass weighted mean, respectively. 



ulation, the total flux of 850 /im in the LBG phase is 
5*850 = 0.65 ~ 0.3 mjy for a galaxy at z = 3. Some 
lensed LBGs have been detected at sub- millimeter, and show 
5g50 = 0.40 mJy and SFR ~ 10 Mq yr"^ at 2 = 2.9 
(MS0451-a : iBorvs et all l2004l). gs.^in = 0.3 9 mJy and 
SFR ~ 24 Mq yr-i (cB58 : iBaker et al.ll2004l ). or Ssso = 
0.65 mJy and SFR ~ 9 M© yr"^ at z = 2.5 (SMMJ16359 
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Figure 8. Total flux at 850 fim in the observed frame as a func- 
tion of evolution time. We set the model galaxy at each redshift 
(circle: z = 2, triangle: z = 3, square: 2 = 4, cross: z = 5 and 
star: z = 6). 
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: iKneib et alJ 1200511. Our results match these observations 
well. Also, Stanwav et al.l (|2010l ) showed the upper limit of 
sub-millimeter flux, 5870 < 0.85 mjy, for LBGs at z — 5. 
Our prediction is under their upper limit, and therefore 
higher sensitivity is required to detect sub-mm flux from 
LBGs at z > 5. On the other hand, there are a few bright 
LBGs detected at sub-millimeter wh ich show 5*850 --^ 1 — 5 
mJy and SFR ^ 100-200 Mp yr~^ jB clokur ov et al.|[2007l : 
ICoppin et"all l2007l : IChapman fc Casev i i2009l' ). These may 
be more massive and dusty star-burst systems. 

In Fig. [9l the map of the brightness temperature at 
A = 850 /im is shown for the model galaxy placed at z = 3. 
The brightness temperature Tb is given by 



Tb = 11.2 



V 850^01 



1" X 1" 



(9) 



The higher Tb traces the star-forming regions. Clumpy 
clouds with higher Tb are distributed over an extended area 
in the LAE phase (tagc ~ 0.1 — 0.3 Gyr), while those are con- 
centrated near the center in the LBG phase (tagc ~ 0.5 — 1.0 
Gyr). 

The upper panels in Figure [9] show Tb in the original 
resolution of simulations that corresponds to ~ 0.13". The 
model galaxy exhibits Tb ~ 0.1 — 40 mK, where the high Tb 
regions are of > 20 mK. The integration time required to 
detect extended sources of ~ 20 mK with the angular res- 
olution of ~ 0.13" is assessed to be ~ 5 hours per each 
beam. The sub-millimeter flux strongly correlates with the 
stellar distribution. At an early phase, the area to cover the 
half mass of stars is ~ 13 arcsec^ . At a later phase, the model 
galaxy becomes as compact as ~ 4 arcsec^ . Hence, the num- 
ber of beams to cover the area will be quite large. Even for 
the later compact phase, the required time to cover the area 
will be a few hundred hours. Thus, the observations with 
the angular resolution of ~ 0.13" do not seem achievable to 
trace the sub-mm structure for early star-forming galaxies. 

The lower panels of Figure [5] show Tb in the coarse- 
grained resolution of 9 ^ 0.65". The high Tb regions show 
> 3 mK. In such an angular resolution, the sensitivity with 
16 antennas of ALMA is ~ 2 mK by ~ 60 minute integra- 
tion. If an early star-forming galaxy is observed, some beams 
may detect sub-mm flux by the integration of less than 1 
hour. Thus, the detection of clumpy sub-mm features of an 
early star-forming galaxy seems feasible. Moreover, obser- 
vations with lower angular resolution of > 1" allow us to 
study the statistics of high-z sub-mm sources, e.g., the lu- 
minosity function of sub-mm galaxies. ALMA will detect 
the sub- millimeter flux from LBGs at 2 ~ 3 by ~ 10 minute 
integration with 16 antennas. Even for a galaxy at 2 = 6, 
ALMA can detect it by ~ 60 minute integration. 

Recently, Dayal et al. (2011) have estimated a bit lower 
flux at 850 fim than the present prediction. They have as- 
sumed homogeneous distributions of stars and dust for an 
LAE identified in a cosmological simulation. Compared to 
our results, it implies that the dumpiness of star-forming 
regions and interstellar medium significantly contributes to 
the enhancement of the sub-mm brightness. 



5 SUMMARY 

We have performed three-dimensional radiative transfer cal- 
culations on high-resolution hydrodynamic simulations with 
inhomogeneous metal enrichment. Then, we have explored 
the dust temperature and sub-mm flux in a high-redshift 
star-forming galaxy. Attention has been concentrated on the 
sub-mm properties of LAE and LBG phases. The star for- 
mation rate (SFR) is ~ 10 M0 yr~^ in the LAE phase, and 
several M0 yr~^ in the LBG phase. As a result, we have 
found that dust grains concentrated on clumpy star-forming 
regions can effectively absorb UV radiation from stars in the 
LAE phase, and then the grains are heated up to Tduat ^ 35 
K. On the other hand, in the LBG phase, the continuous 
supernovae blow away dusty gas from star-forming regions. 
Hence, the grains cannot absorb much radiation from stars 
and turn into a cold state. Resultantly, the emergent sub- 
millimeter brightness is peaked in an LAE phase around 
fagc ~ 0.3 Gyr, independent of the assumed redshift. This 
shows that the sub-millimeter brightness does not neces- 
sarily increases according to the increase of the dust-to-gas 
ratio or the metallicity. We have found that the dumpiness 
of star-forming regions and interstellar medium significantly 
enhances the sub-mm brightness. 

Also, by deriving the spectral energy distributions, we 
have assessed the flux at 850 /im. The flux is found to be 
5850 ~ 0.2 — 0.9 mJy, if the model galaxy is placed at 

6 ^ 2 ^ 2. Even for 2 = 6, the sub-mm flux does not 
decrease largely because of negative K-correction. With an- 
gular resolution of 9 ^ 0.65", ALMA can detect the sub-mm 
flux from an LAE as well as an LBG by the integration of 
less than 1 hour with 16 antennas. Observations with angu- 
lar resolution of > 1" allow us to detect high-z sub-mm 
sources by ~ 10 minute integration with 16 antennas of 
ALMA. Therefore, we can study clumpy sub-mm features 
of an early star-forming galaxy and also the statistics of the 
luminosity function of sub-mm galaxies with ALMA. 
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